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ABSTRACT 

In the absence of H2 molecules, the primordial gas in early dark matter halos with 
virial temperatures just above Tvir ^lO^K cools by coUisional excitation of atomic H. 
Although it cools efficiently, this gas remains relatively hot, at a temperature near T ^ 
8000 K, and consequently might be able to avoid fragmentation and collapse directly 
into a supermassive black hole (SMBH). In order for Il2-formation and cooling to be 
strongly suppressed, the gas must be irradiated by a sufficiently intense ultraviolet 
(UV) flux. We performed a suite of three-dimensional hydrodynamical adaptive mesh 
refinement (AMR) simulations of gas collapse in three different protogalactic halos 
with Tvir >IO'*K, irradiated by a UV flux with various intensities and spectra. We 

determined the critical specific intensity, J™*, required to suppress H2 cooling in each 
of the three halos. For a hard spectrum representative of metal-free stars, we find (in 
units of 10"^^ erg s~^ Hz"^ sr~^ cm~^) 10'* < J|"* < 10^, while for a softer spectrum, 
which is characteristic of a normal stellar population, and for which H~ -dissociation 
is important, we find 30 < J2"' < 300 . These values are a factor of 3-10 lower 
than previous estimates. We attribute the difference to the higher, more accurate H2 
collisional dissociation rate we adopted. The reduction in J|"* exponentially increases 
the number of rare halos exposed to super-critical radiation. When H2 cooling is 
suppressed, gas collapse starts with a delay, but it ultimately proceeds more rapidly. 
The infall velocity is near the increased sound speed, and an object as massive as 
M ^ 10^ Mq may form at the center of these halos, compared to the M ~ 10^ Mq 
stars forming when H2 -cooling is efficient. 

Key viTords: cosmology:theory - black holes physics - methods:numerical 



1 INTRODUCTION 

The discovery of very bright quasars, with luminosities 
> 10*^ erg s~*, at z ~ 6 in the Sloan Digital Sky Survey 
(SDSS) suggests that some SMBHs as massive as a few times 
10^ Mq already existed when the universe was less than 1 
Gyr old (see, e.g.. Fan 2006 for a review). In principle, these 
large black hole (BH) masses, inferred from the apparent 
luminosities, could have been overestimated due to strong 
gravitational lensing and/or beaming. However, no obvious 
sign of either effect was found in the images or spectra of 
these quasars (Willott et al. 2003; Richards et al. 2004). 
Explaining how SMBHs with masses in excess of 
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10^ Mq could assemble within 1 Gyr presents some chal- 
lenges. Perhaps the most natural proposal is that they grow, 
by a combination of Eddington-limited accretion and merg- 
ers, from the stellar-mass seed BHs provided by the rem- 
nants of the first generation of massive, metal-free stars (e.g., 
Haiman & Loeb 2001). Indeed, the initial seed BHs, with 
masses of the order of their progenitor stars, ~ 100 Mq 
(Abel et al. 2000, 2002; Bromm et al. 2002), are expected to 
be present very early (at redshifts prior to z >10). There 
are, however, a number of potential difficulties with this 
scenario. First, the early seeds must accrete near the Ed- 
dington rate for a Hubble time, without any prolonged in- 
terruption (Haiman & Loeb 2001), which requires that the 
BHs are continuously surrounded by dense gas (e.g. Turner 
1991; Alvarez et al 2009). However, early seed BHs are ex- 
pected to undergo frequent mergers, and the gravitational 
waves emitted during the BH mergers impart a strong re- 
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coil to the coalesced BHs (e.g. Pretorius 2005; Campanelli 
et al. 2006; Baker et al. 2006). The typical velocity of this 
gravitational kick is expected to exceed ~ 100 km/s, which 
is significantly higher than the escape velocity ( ;^10 km/s) 
from typical dark matter halos at 2 ~ 10. BHs are therefore 
easily ejected, or at least displaced from the dense nuclei 
of their host halos at high redshift, interrupting their ac- 
cretion (Haiman 2004; Yoo & Miralda-Escude 2004; Shapiro 
2005, Volonteri & Rees 2006; Blecha & Loeb 2008; Tanaka & 
Haiman 2009). Even if such disruptive kicks are avoided (be- 
cause mergers at early times may be rare and/or occur pri- 
marily between unequal-mass BHs; Volonteri & Rees 2006; 
Tanaka & Haiman 2009), there remain two fundamental dif- 
ficulties. First, when the effects of photoheating and radi- 
ation pressure are included, the maximum allowed steady 
accretion rate is significantly reduced, at least in spherical 
symmetry (e.g. Milosavljevic et al. 2009a), suggesting that 
accretion must be intermittent, with a time-averaged rate 
well below the Eddington-limit (e.g. Ciotti & Ostriker 2001, 
2007; Milosavljevic et al. 2009b). Second, models in which 
sufficient numbers of early BHs are able to accrete and grow 
to 10^ M0 by 2: « 6 tend to overproduce the abundance of 
« 10^ Mq BHs by several orders of magnitude (compared 
to the abundance inferred from local observations). This re- 
quires a preferential suppression of BH growth in low-mass 
halos, plausibly due to negative feedback effects acting on 
these halos (Bromley et al. 2004; Tanaka & Haiman 2009). 

An alternative way of assembling SMBHs is through 
more rapid (super-Eddington) accretion or collapse. In this 
family of models, primordial gas collapses directly into a BH 
as massive as lO" - 10*^ M© (Oh & Haiman 2002 [hereafter 
OH02]; Bromm & Loeb 2003 [hereafter BL03]; Koushiap- 
pas et al. 2004; Lodato & Natarajan 2006; Spaans & Silk 
2006; Begelman et al. 2006; Volonteri et al. 2008), possi- 
bly onto a pre-existing smaller seed BH (Volonteri & Rees 
2005), or through the intermediate state of a very massive 
star (BL03). Such a head-start evades problems encountered 
by the models where SMBHs grow at the Eddington limit 
from stellar mass seeds. A necessary condition for such direct 
collapse models is that the collapsing gas avoid fragmenta- 
tion; otherwise, normal Pop HI stars would be producedQ 
At the density of ~ lO'^cm""^ (the critical density for H2; see 
below), the Jeans mass is Mj ^ 10*^ Mq {T/10* K)^/^ The 
gas temperature T depends crucially on whether H2 cooling 
is efficient: T « 100 K, achievable if H2 cooling is efficient, 
would imply Mj « 10'^ Mq , so that PopIII stars might form, 
whereas T « 10* K, expected in the absence of H2, would 
yield Mj ^ 10^ M©, suggesting that direct collapse into a 
M « 10" Mo SMBH may be feasible (OH02). 

Numerical simulations have indeed shown that fragmen- 
tation is inefficient when H2 cooling is absent (BL03; Regan 
& Haehnelt 2009a, 2009b). However, in most models, the ab- 
sence of H2 was assumed, rather than derived. The notable 
exceptions are Spaans & Silk (2006), whose model does not 
require any explicit H2 destructiorjj, and BL03, who per- 

^ Another necessary condition is for the gas to loose angular mo- 
mentum efficiently. Other than the angular momentum transfer 
occurring above the resolution of our simulations, this topic will 
not be addressed in the present paper. See, e.g. Begelman et al. 
(2006) for a discussion and for references. 

^ Spaans & Silk propose that over a relatively narrow range of 



formed simulations with an H2-photodissociating Lyman- 
Werner background. The absence of H2 molecules from pro- 
togalactic halo gas can be justified by a sufficiently intense 
UV radiation, either in the Lyman-Werner bands, directly 
photo-dissociating H2 (near a photon energy of ~ 12eV) 
or photo-dissociating the intermediary H~ (photon energies 
^0.76eV). The relevant criterion is that the photodissocia- 
tion timescale is shorter than the H2-formation timescale; 
since generically, tdiss oc J and tform oc p, the condition 
tdiss = tform yiclds a critical flux that increases linearly with 
density, J'^'^'* oc p. In "minihalos", with virial temperatures 
Tvir < lO^K, the gas cannot cool in the absence of H2, the 
densities remain low (~ lcm~^; e.g. Mesinger et al. 2006) 
and H2 can be dissociated even by a relatively feeble UV 
flux. The critical value has been found to be J21 ~ 0.1 
(Haiman et al. 1997; Machacek et al. 2001, 2003; Mesinger 
et al. 2006, 2008; Wise & Abel 2007; O'Shea & Norman 
2008; here and in the rest of the paper, J21 denotes the 
speciflc intensity just below 13.6eV, in the usual units of 
10~^^erg cm~^ sr~^ s~^ Hz~^). This value is much smaller 
than the expected level of the cosmic UV background in the 
Lyman-Werner bands near reionization (BL03) 

1 he N-^YuPh 

Jbg ~ -;; — (^) 

/esc 47r rrip 

or 

where /osc is the escape fraction of ionizing radiation, A^^ is 
the average number of photons needed to ionize a hydrogen 
atom, Yh = 0.76 is the mass fraction of hydrogen, nip is the 
proton mass, and pb is the background baryon density with 
Qth'^ = 0.023. 

The critical intensity J21'' in larger halos, with virial 
temperatures T^ii ^ 10* K, however, is much higher 
(Omukai 2001, hereafter OMOl; OH02; BL03). This is pri- 
marily because the gas in these halos can cool via excitations 
of atomic H and reach much higher densities, and because 
the H2 molecules can then become self-shielding (OH02). In 
particular, for halos with Tvir ~ 10*K, the value has been es- 
timated in one-zone models to be J^i^ « 10^ - lO'^ (OMOl). 
This range covers different assumed spectral shapes; in par- 
ticular, a thermal spectrum with T. = 10* - lO'^K (OMOl). 
Using three-dimensional smooth particle hydrodynamics 
(SPH) simulations, BL03 find J2T' ^lO'^ (for T, = lO^K) 
and Jaf* ^10^ (for T. = 10*K), in agreement with the one- 
zone results. 

The UV background that illuminates collapsing Tvir ~ 
10*K halos was likely established by the massive pop-Ill 
stars that had formed in previous generations of miniha- 
los. In this case, the background spectrum is likely to be 
closer to the T« = lO^K case, implying J^i*" gsimlO^ . Fur- 
thermore, the early minihalos are expected to be easily 
self-ionized, with most of their ionizing radiation escap- 
ing into the intergalactic medium, i.e. /esc ~ 1 (Kitayama 
et al. 2004; Whalen, Abel & Norman 2004). Equation [2] 

densities and hydrogen column densities, the hya photons emit- 
ted by atomic H cooling are trapped within the collapsing gas - 
this prevents the temperature from falling below ~ 8, OOOK, and 
keeps the H2 molecules coUisionally dissociated. 
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shows that, unless /esc is much smaUer, and/or A'^. is large 
iUc/N^ ^10-^), the mean cosmic background is unlikely 
to reach the required critical value. The background will 
have inevitable spatial fluctuations, and a very small fraction 
/(J > J"^^) of 10'' K halos that have an unusually close and 
bright neighbor may still see a sufflciently high flux. Dijkstra 
et al. (2008) used a model for the three-dimensional spa- 
tial clustering of halos to estimate this fraction, and found 
/(J > 10^) ~ 10~^, with an exponential dependence of this 
result on Jal^EI 

In this paper, we derive detailed estimates for J21'* 
based on a suite of 3D hydrodynamic simulations. Our moti- 
vation is two-fold. First a UV intensity that exceeds J21'' is 
crucial for the feasibility of direct SMBH formation models. 
Second, the existing estimates of J21'' significantly exceed 
the expected value of the mean cosmic background. If these 
estimates are correct, then J > J21'* will be experienced only 
by those rare halos that probe the bright tail of the spatially 
fluctuating background J. In this case, even a small change 
in the value of JJ"' can cause a large change in the expected 
number of halos that can form SMBHs by direct collapse. 

Our paper adds to the earlier work of OMOl, which 
estimated J21'* based on a one-zone model with a fixed, 
prescribed collapse dynamics, that could not address frag- 
mentation, and to the results of BL03, who use three- 
dimensional simulations, but report only approximate upper 
and lower limits on J^f* for a single halo. In addition, we 
study the behavior of the collapsing gas in detail as a func- 
tion of J21 , and obtain a rough estimate for the fluctuations 
in J|"' by following the collapse of three different halos. We 
also obtain a quantitative estimate of the final collapsed cen- 
tral massive object, based on the infall time-scales observed 
in each case. 

The rest of this paper is organized as follows. In § [51 we 
describe the simulation setup in detail. In § [S] we present 
the range of J21'' values derived from the simulations, and 
explain the underlying physics using a one-zone model. We 
also estimate the fluctuations in J21'' . In § [l] we discuss the 
ultimate fate of the halos, for different values of J21. Finally, 
we summarize the results and present the conclusion of this 
work in § [5] 



2 NUMERICAL METHODOLOGY 

We use the Eulerian adaptive mesh refinement (AMR) code 
Enzo, which has been tested extensively and is publicly 
available (Bryan 1999; Norman & Bryan 1999; O'Shea et al. 
2004). Enzo uses an N-body particle-mesh solver to follow 
dark matter dynamics, and an Eulerian AMR method by 
Berger & Colella (1989) to solve the hydrodynamical equa- 
tions for an ideal gas. This combination allows for high dy- 
namic range in gravitational physics and hydrodynamics. 
Nested grids are used whenever higher resolution is needed. 



^ Alternatively, the critical value could be established by sources 
internal to the halo, e.g. by a vigorous phase of starburst (Omukai 
& Yoshii 2003) or by an accreting stellar seed BH (Volonteri 
& Rees 2005). However, having gone through star-formation al- 
ready, the halo gas is unlikely to still be metal— free, and is then 
likely to fragment into low— mass stars, rather than collapsing di- 
rectly into a SMBH (Omukai et al. 2008). 



At each new refinement level, a parent grid is replaced by a 
few smaller child grids. 

The chemical composition of the gas is followed by solv- 
ing the non-equilibrium evolution of nine species: H, H^, 
He, He+, He++, H", Hj, H2, and e" (Abel et al. 1997; 
Abel et al. 2000). Our reaction network did not include HD 
or other molecules involving deuterium. This should make 
very little difference to our results, since HD cooling only 
becomes important for temperatures below a few hundred 
Kelvin (McGreer & Bryan 2008). The H2 radiative cooling 
function of Galli & Palla (1998) is employed to follow the 
temperature of the gas. 

A few modifications were made to the publicly available 
Enzo code. First, we added direct H~ photodissociation into 
chemistry solver. 



H' 



hi^ H-l-e 



(3) 



As explained below, this reaction is important in determin- 
ing the value of J21'' . Second, we included the self-shielding 
of H2 in the LW bands when computing H2 photodissocia- 
tion rate. Specifically, the intensity in the LW band is multi- 
plied by a self-shielding factor fsh given by Draine & Bertoldi 
(1996), 



fsh 



Ho 



10"cm- 



-3/4 



(4) 



where A^H2 is the H2 column density. Since A'h2 is a non- 
local quantity, it is computationally very expensive to obtain 
its exact value. To save computing time, we made the com- 
monly used approximation. 



/H2'^totAj, 



(5) 



where , "-tot and Aj are the H2 fraction by number, the 
total particle number density, and the Jeans length, respec- 
tively. With this approximation, and fsh can be com- 
puted from the local values of the temperature, density and 
H2 fraction. We will quantify the accuracy of this approxi- 
mation in i) 13.51 below. 

The simulation is set up in a comoving box of size 
1 Mpc, assuming a standard ACDM model with the fol- 
lowing parameter: f^DM = 0.233, fit, = 0.0462, ~ 0.721, 
as = 0.817, Us = 0.96 and h = 0.701 (Komatsu et al. 2009). 
We first perform a preliminary run in order to identify halos 
suitable for detailed study. This run has a root grid with a 
resolution of 128"^ and no nested grids. Radiative cooling is 
turned off, so that the gas in the halos is unable to contract 
to high densities in this run. We evolved the simulation to 
z = 10, stopped it, and used the HOP halo finding algorithm 
(Eisenstein & Hut 1998) to identify dark matter halos in the 
output files. Three halos, which are labeled below as A, B 
and C, were chosen for high-resolution re-runs. All three of 
these were selected to have virial masses of a few X 10^ Mo 
at a redshift of z; = 10. In this paper, virial mass is defined 
to be the total mass, including both gas and dark matter 
components, inside a spherical averaged overdensity of 200 
with respect to the critical density of the universe. 

We generated a new set of initial conditions for the three 
chosen halos. Three nested grids with twice finer resolution 
were added, so that the effective resolution of the innermost 
grid was 1024^, resulting in a dark matter particle mass of 
86 Mq. Radiative cooling was turned on in the re-runs, and 
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the grid cells were adaptively refined based on the follow- 
ing three criteria: baryon mass, dark matter mass and Jeans 
length. According to the first two criteria, additional grids 
are added when the baryon (dark matter) mass in a grid cell 
exceeds 68 (683) solar mass, corresponding to 4 (8) times the 
initial mass in one grid cell (particle) in the most refined re- 
gion. The third criterion ensures that the Jeans length is 
resolved by at least 4 grid cells, so that no artificial frag- 
mentation would take place. In addition, to avoid numerical 
effects due to the finite mass of dark matter particles, the 
gravity of dark matter particles is smoothed at refinement 
level 13, which corresponds to a smoothing scale of 0.954/^ 
(comoving) parsec. Each dark matter particle has a mass of 
~ 85 Mq . We allow the simulations to proceed until a max- 
imum refinement level of 18 is achieved, corresponding to a 
resolution of 0.0298//i (comoving) parsec, or about 800 AU 
(absolute). 

For each of the three halos, we ran a series of simulations 
with different UV spectra and intensities. For the spectral 
shape, we adopted a Planck spectrum with a black-body 
temperature of either T* = 10* K or T. = lO'^ K (hereafter 
denoted by T4 and T5, respectively). The softer of these 
spectra is meant to approximate the mean spectrum of a 
normal stellar population, whereas the higher-temperature 
case is closer to the harder spectrum expected to be emitted 
by the first generation of massive, metal-free stars (Tumlin- 
son & ShuU 2000; Bromm, Kudritzki & Loeb 2001; Schaerer 
2002). Using these two spectral types allows us to compare 
our results with previous work (OMOl, BL03) which adopted 
the same spectral shapes. In Table [1] ([2]), we list the red- 
shift (^coi), virial mass (rrivir.coi) and central gas tempera- 
ture (Tcont) in the halos when their cores collapse in the 
presence of type T4 (T5) UV background. Here, "core col- 
lapse" is simply defined as the time when the maximum 
refinement level (level 18) is reached. In practice, once the 
collapse starts, it proceeds very rapidly. As a result, the 
refinement level adopted for this definition makes little dif- 
ference to our results, as long as it is chosen to be at level 
13 or higher. We generally varied J21 by factors of 10, but 
included additional runs with J21 = 3 x 10^ for halos B and 
C in the T4 case in order to determine Jfi'' more precisely. 
The values for the J21 = 10'' case for halo C are missing 
from Table [TJ because the halo moved out of the refinement 
region before it collapsed in this run. This occurred because 
of the late collapse redshift for this halo (recall that halos 
were selected a.t z — 10). We could have rerun the simula- 
tion with a larger refined region, but this was unnecessary 
because the J21 = 3 x 10^ run was sufficient for determining 

Jcrit 
21 ■ 



3 THE CRITICAL VALUE OF J21 

3.1 Results from Three— Dimensional Simulations 

Based on the central gas temperatures (Tcont) listed in Ta- 
bles [T] and [21 the halos can be grouped unambiguously into 
two categories: "cool" halos and "hot" halos. The "cool" 
halos all have central temperatures near ~ 1000 K, whereas 
the "hot" halos all have central temperatures near ~ 6000 
K. To see the differences between these two categories more 
clearly, we show in Figures [T] - [4] the spherically averaged 
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Figure 1. Spherically averaged profiles of the particle density, gas 
temperature, e~ fraction and H2 fraction in halo A, for different 
values of the intensity J21 of a type T4 background. 
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Figure 2. Same as Figure [T] except for type T5 backgrounds. 
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Figure 3. Same as Figure [T] except for Halo C. 



Keeping Protogalactic Gas H2 Free 5 



Table 1. Redshift (zcol)i virial mass (mvir,col) and central gas temperature (Tccnt) of the halos when their cores collapse in the presence 
of a T* = 10-*K black-body ("T4") UV background. 
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Table 2. Same as Table[T] but for a T, = lO^K black-body {"T5") UV background. 
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the "cool" halos start to cool via H2, 
stay hot. As the figures show, R ~ 



while the "hot" halos 
10 pc corresponds to 
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Figure 4. Same as Figure |3] except for type T5 backgrounds. 



profiles of the particle density, gas temperature, e~ fraction 
and H2 fraction in halos A and C, for different values of J21, 
and in the presence of type T4 and T5 UV backgrounds. The 
profiles of halo B are not shown because they are very sim- 
ilar to those of halo C. All profiles are shown at the time of 
collapse at ^coii- The collapse is delayed significantly as the 
intensity of the background is increased (see Tables [1] and 
so that the curves shown in each panel reflect conditions 
at significantly different redshifts. 

Consistent with the values in Tables [1] and [2j these fig- 
ures show that central temperatures are very different for 
"hot" and "cool" halos. However, the temperature profiles at 
large radii are actually similar in the cases where J21 ^ 10. 
As the gas falls in, the temperature first rises due to viri- 
alization shocks, and stays at 8000 K because of atomic 
cooling. The clear differences arise around i? 10 pc, where 



a density of n ~ 10^ cm"''. The conditions at this density 
are crucial in determining the value of J21''. In particular, 
the difference in the temperature profiles is clearly induced 
by the difference in the H2 fractions. As the bottom right 
panel of each of Figures [T] - |4] shows, the H2 fractions in- 
crease with decreasing radii in the "cool" halos. In the "hot" 
halos, the fractions increase with decreasing radii outside 
J? 10 pc, but the fraction stays constant, at the low level 
of ~ 10~*, inside this radius. There are also clear differences 
in e~ fraction profiles: the cores of the "hot" halos have a 
much larger electron fraction (~ 10~^), whereas in the cores 
of the "cool" halos, the electrons are depleted. Finally, al- 
though the density profiles are overall quite similar, the gas 
densities reached within the inner 10 pc in the "hot" halos 
are noticeably higher than in the "cool" halos. 

As we will discuss in more detail in § U below, dur- 
ing the process of quasi-static contraction, the density and 
temperature adjust so that the cooling, sound crossing, and 
dynamical time scales are all approximately equal (tcooi ~ 
Across ~ iff ). From the condition tcross ~ iff , we find nR^ oc T 
(see equation 1161 below), which explains why the density is 
higher when temperature is higher at a given radius. Within 
J? ~ 10 pc, the e~ fractions in "hot" halos have a much shal- 
lower slope than in "cool" halos. This is because the main 
reactions determining the electron fraction (reactions 1 and 
4 in the Appendix) are highly temperature sensitive - in par- 
ticular, recombinations (reaction 4) are exponentially more 
rapid at lower temperatures, and the gas is therefore much 
more fully recombined in the "cool" halos. 

It is worth noting that in the case of J21 — 1 and type 
T5 UV background, the halo gas is able to collapse before 
the halo potential grows deep enough to reach a virial tem- 
perature of 8000 K, where atomic cooling becomes im- 



6 C. Shang, G. L. Bryan and Z. Haiman 



Table 3. Jg'* for the three different halos (A, B, and C), and 
for the two different UV background spectra (T4 and T5). 





Halo A 


Halo B 


Halo C 


One-zone 


T4 


3 X IQi - 10^ 


10^ - 3 X 10^ 


10^ - 3 X 10^ 


3.9 X IQl 


T5 


10"' - 3 X 10"' 


3 X 10" - 10^ 


3 X 10"* - 10^ 


1.2 X 10* 



portant. Efltectively, at this low UV background, the halos 
are "minihalos" , whose properties are fully determined by 
H2 cooling. This causes a "break" in the properties of halos, 
such as their mass and collapse redshift, between J21 = 1 
and J21 = 10. Such a "break" was also seen in O'Shea & 
Norman (2008), but at J21 ~ 0.1. The fact that this break 
occurs at a higher J21 in our simulations is not surprising, 
since we include self-shielding, which was ignored in O'Shea 
& Norman (2008). 

Finally, we note that even in the cases with higher 
values of J21, when the collapse is delayed, and the dark 
matter halo potentials grow somewhat deeper, the effective 
virial temperature of the halos we examined are still only 
marginally (if at all) above IC^K. Indeed, Figures [T]-|4] show 
that the gas in the simulated halos is never fully collision- 
ally ionized - the highest free electron fraction reached is 
only 10""^. This regime differs from the common assumption 
that "second generation" halos are coUisionally ionized (e.g. 
OH02). In the strict absence of any H2 cooling, the shocks 
that occur in more massive halos could, of course, eventu- 
ally produce full ionization. However, prior to building up 
this higher mass, every halo must go through the "marginal" 
stage where its effective virial temperature is close to, but 
just above, Tvir ~ 10'*K. As our results show, the gas at this 
stage can already cool efficiently via excitations of atomic 
H. Therefore, we expect that halos that reach virial temper- 
atures significantly above « lO^'K and become fully ionized, 
but still consist of pure atomic H (and He), will be exceed- 
ingly rare. 

Because each of our runs produce a collapsed halo that 
belongs unambiguously either to the "cool" or "hot" cate- 
gory, we simply define J21'' as the value that divides these 
two regimes. Our suite of simulations is then sufficient to 
determine 72"* to within a factor ~ 3 for each of the three 
halos, and for the two types of UV spectra. The results are 
listed in Table [3] 

There are two interesting new points that can be con- 
cluded from the J21'' values in Table [S] First, the ranges of 
J2"* we found are smaller than previously estimated. Pre- 
vious studies (OMOl; Omukai et al. 2008) found J^j'* 10^ 
for the type T4 spectrum and ~ 3 x 10^ for the type T5 
spectrum. The upper and lower limits of ^10^ and ^10^, 
reported from SPH simulations by BL03 for these two types 
of spectra, are consistent with the above values. These val- 
ues, however, are a factor of ~ 10 larger than our results for 
halo A, and a factor of ~ 3 larger than for halos B and C. 
Second, even with our crude sampling of J values, we see a 
"scatter" in J21'': the critical flux is a factor of '--^ 3 lower for 
halo A than for halos B and C, despite the similar masses 
and collapse redshifts of these halos. 

The fact that we find J21'' values that are smaller than 



previously estimated is particularly important, because the 
critical UV fluxes are high compared to the expected level 
of the cosmic UV background at high redshifts. As noted 
above, however, the background will inevitably have spatial 
fluctuations, and a small fraction of halos may still see a suf- 
ficiently high flux. Dijkstra et al. (2008) used a model for the 
three-dimensional spatial clustering of halos to estimate the 
probability distribution function (PDF) of the background 
LW flux Jbg, as sampled by DM halos with Tvir « lO^K. 
Their results show that the interesting range of J21'* samples 
the bright, steeply falling tail of the flux PDF. In particular, 
in their Figure 2, Dijkstra et al. (2008), show the fraction 
/(> Jl"*) of halos exposed to a flux above a given J|"* . This 
fraction is very sensitive to J^f': for J^f' ^lO"*, the PDF 
drops to negligibly low values (/ ^10"**); for Jg'* = 10^, 



/ 



10" 



whereas for Jg'* = 10% / ~ 10"^ This im- 



plies that the 10-fold decrease we find in the value of JJ"* 
increases the number of candidate DM halos, where direct 
SMBH formation may be feasible, by a factor of ~ 10^. 

3.2 Results from One— Zone Models 

In order to understand the above results, we performed one- 
zone calculations similar to those in OMOl. The chemical 
reaction and cooling rates are set to be the same as in the 
Enzo runs, but the gas density is assumed to have a single 
value that follows a fixed prescribed evolution. The list of 
chemical reactions we used, and their rates, are shown in the 
Appendix. In this model, the dark matter and the baryons 
start their collapse at the turnaround redshift, which is set 
to be 2; = 17. At this redshift, zero initial velocities are 
assumed. The subsequent evolution of the dark matter den- 
sity is computed using the spherical collapse model, up to 
the time of virialization, after which it is assumed to stay 
constant at the virial density. The evolution of the baryonic 
component is followed using the equation, 

dp;, _ dpi, 
dt ^ ts " 



(6) 



where pb is the baryonic density, and tg — y%37r732Gp is 
the dynamical (free-fall) time. The thermal evolution is de- 
scribed the following equation, 

de _ _ _d_ / J_\ _ Anct 
dt ^ ^dt \pbj Pb ' 

where e is the internal energy per unit baryonic mass, p is 
the gas pressure, and Anct is the net cooling rate computed 
using cooling function by Galli and Palla (1998). The inter- 
nal energy density is computed from 

1 kbT 



7 - 1 p.mH ' 



(8) 



where 7 = 5/3 is the adiabatic index, fee is the Boltzmann 
constant, p is the mean molecular weight, and mn is the 
mass of a hydrogen nucleus. We include H2 self-shielding in 
the same way as described above for the Enzo runs. 

Figure [5] shows the gas temperature, electron fraction 
and H2 fraction as a function of density computed for halo 
A with Enzo (thick curves) and in the one-zone model (thin 
curves) for both type of UV backgrounds and for both "cool" 
halos (solid curves) and "hot" halos (dashed curves), with 
J21 values as indicated in the middle panels. 
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Figure 5. Gas temperature, electron fraction and H2 fraction as a 
function of density computed for halo A witli Enzo (tliick curves) 
and in a one-zone model (thin curves). The solid and dashed 
curves show result for "cool" and "hot" halos, respectively, with 
the values of J21 as indicated in the middle panels. The left panels 
are for type T4 UV backgrounds, and the right panels are for type 
T5 UV backgrounds. 



The results from the simulation and the one-zone model 
are in excellent agreement at n 10^ cm~^. There is a clear 
difference in the temperatures at n ^ 10^ cm~^, which is 
expected, since the one-zone model prescribes a smooth adi- 
abatic collapse, and the heating of the gas is purely due to 
this adiabatic compression. In contrast, in the simulation, 
the gas experiences shocks, which elevate the temperatures 
in the low-density regime. To check whether the shocks af- 
fect the comparison at higher densities, we mimicked the 
shocks in the one-zone model by artificially setting the tem- 
perature to ~ 8000 K at a fixed low density. We found that 
this makes the evolution match the simulations better at 
n < lO'^ cm~^, but has little effect at higher density. In par- 
ticular, JJf' computed in the one-zone model (see below) 
changes by ^ 3%. This is because, as already mentioned 
above (and will be discussed further below), J21'' is primar- 
ily determined by the conditions at the critical density of 
H2, which is ~ 10^ cm"^. 

We performed one-zone calculations similar to those 
shown in Figure (5] but varied J according to a Newton- 
Raphson scheme, until we converged on the critical value. 
We found Jg" = 39 for type T4 backgrounds, and J|f = 
1.2 X 10* for type T5 backgrounds. These values fall within 
the range of J21'' identified in the Enzo runs for halo A, but 
are slightly below the Enzo range for halos B and C. We will 
discuss possibly sources of this difference below. 



10"' 



10"' 



10-^' 



10"' 



1 1 1 1 1 III 

\ 


1 i i 1 i 


1 1 1 i i 1 i II 

Sfgrm 
fed 


N 

\ 

\ 

\ 

\ 

\ 

\ N 
\ \ 
\ \ 












0^ yX 





10" 



10' 



10' 



lO'' 
n(cm"'') 



10' 



10= 




-'form 
Sod 
cfpd 
-^totd 



10' 



lo-' 



n(cm"'') 



10° 



Figure 6. The normalized reaction speed 5 as a function of n 
for type T4 backgrounds. Solid curves: formation speed Sform- 
Dotted curves: photo-dissociation speed Sp^. Dashed curves: col- 
lisional dissociation speed Scd- Dash-dotted curves: total dissoci- 
ation speed Stotd- The thin and thick curves are for J21 = 10^ 
and J21 = 10^, respectively. 



10° 



Figure 7. 5 as a function of n, as in Figure \E\ but for type 
T5 backgrounds. Line types are the same as in Figure |6l except 
here the thin and thick curves are for J21 = 10"* and J21 = 10^, 
respectively. 



3.3 The importance of collisional dissociation 

The one-zone model described above enables us to study 
the effect of each reaction on J21'', thereby understanding 
the physics that determines Jli'': and identifying the uncer- 
tainty in J21'' due to uncertainties in the chemical reaction 
rates. 

The five most important reactions in determining J21'* 
are (OMOl): H2 formation (reaction 10), H~ formation (re- 
action 9), H~ photo-dissociation (reaction 25), H2 photo- 
dissociation (reaction 28) and H2 collisional dissociation (re- 
action 15). Reactions 9 & 10 are the main channel of forming 
H2 molecules. Reaction 25 competes with reaction 10 for H~, 
modulating the effective H2 formation rate (OMOl), 



^form — ^9 



fcion -I- k2 



(9) 
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where fci is the rate of reaction i (as listed in the Appendix). 
Both photo-dissociation and colhsional dissociation compete 
with this effective H2 formation process. The H2 fraction in 
the gas is given approximately by the equilibrium value for 
whichever dissociation process dominates, 



// "'form J. 
H2 = mm /en, 



k2 



fel5 



(10) 



To illustrate how these reactions determine the critical value 
of J21, in Figures |6] and [7] we show the normalized reaction 
speed S as a function of n for different values of J21. The 
normalized formation, photo-dissociation, colhsional disso- 
ciation and total dissociation speed are defined as 



'S'form — kfoj-^TLcTlJl / n ~ ^form./c, 
'S'totd ~ Spd + Scd- 



(11) 

(12) 
(13) 
(14) 



In these figures, the temperature, /hj and /c are set to 
8000 K, 10"^ and 5 x 10"^, respectively (with the self- 
shielding factor /sh computed from equations U and O. 
These choices were motivated by Figures IH4I which show 
that never exceeds ~ 10~^ for the "hot" halos, and /c 
is between 10~^ and 10~^ at n ~ 10^. We therefore adopt 
10~^ as the critical H2 fraction, above which the gas could 
cool below ~ 8000 K. 

The solid curves in Figures |6] and [7] show the formation 
speeds Sform, and the other three curves show the photo- 
dissociation, colhsional dissociation, and total dissociation 
speeds, as labeled. Note that at low density, Sform increases 
with density, but at high density, it stays approximately con- 
stant. This behavior follows from equation ([9]): in the low 
density limit, fcion ^ fes, feform ~ fcgfcion/fes cx n; whereas 
in the high density limit, fcion S> A;25, fcform « fcg ~ constant. 
Likewise, the normalized dissociation rate decreases with 
density when it is dominated by Spd, but asymptotes to 
a constant in the high-density limit, where colhsional dis- 
sociation dominates. Figures |6] and [7] show that Scd transits 
from a low value to a high value around the critical density 
ricr ~ 10^ cm"^. 

Most importantly, the formation and dissociation 
speeds are both functions of J. In particular, as J is in- 
creased, the total dissociation rate is increased, but only 
below the density where H2 photo-dissociation dominates 
colhsional dissociation - the colhsional rates are indepen- 
dent of J. Likewise, as J is increased, the effective forma- 
tion rate decreases due to H~ photo-dissociation, but only 
at low densities, where H~ photo-dissociation is more im- 
portant than colhsional H~ dissociation (since, again, the 
colhsional H~ dissociation rate is independent of J). 

The critical value, J"^^ , can be identified as the lowest 
value for which the formation speed and the total dissocia- 
tion speed can become equal, for any value of the density. If 
J is lower than this J"^^ , then, by definition, dissociation is 
always more rapid than formation, and the H2 abundance 
can not increase above the critical value. However, if J ex- 
ceeds this J'^'^'', then at the density where Stotd ~ Storm, the 
H2 abundance increases, and the gas starts to cool. Because 
Scd drops further as the gas cools (see Fig. [S] below), more 
H2 molecules are able to form, and the cooling proceeds in 
a runaway fashion. 



Since the gas initially starts at low density, and, as 
shown in Figures [6] and [T] S'totd is high at both low and 
high densities, whether or not a halo would cool is mainly 
determined by the condition around the density nmin where 
Stotd reaches a minimum. The value at this minimum, 
Wmin ~ 10^ — 10* cm~^ when J21 is close to J2f^, coin- 
cides with the critical density for H2, Tier ~ 10^ cm"'' (see 
below). 

Figures [6] and [7] also show the relative importance of 
H~ photo-dissociation versus direct H2 photo-dissociation 
in setting the critical value of J21. In particular, the J- 
dependence of Sform (solid curves) arises from H~ photo- 
dissociation, whereas the J-dependence of S'totd (dash- 
dotted curves) is from H2 photo-dissociation. Comparing 
the thick and thin solid curves in Figures |6] and [T] we see 
that around the critical density of ~ 10^ cm~^, H~ photo- 
dissociation is much more important in the T4 case. This 
is primarily because of the softer shape of this spectrum: 
for a fixed fiux at 13.6eV, the flux at ~ leV, just above 
the H~ photo-dissociation threshold, is much larger for the 
T4 than for the T5 case (see Omukai et al. 2008 for an ex- 
plicit comparison). Interestingly, in the case of the T4 spec- 
trum, as J is varied, the minimum in Stotd moves almost 
in parallel with the Storm vs. n curve. As a result, J21'' is, 
in fact, quite insensitive to the direct H2 photo-dissociation 
rate (and therefore also to the details of our treatment of 
self-shielding). We verified this conclusion explicitly, by ar- 
tificially setting fc28 = in the one-zone model. In this case, 
J21'* only increases by a factor of ~ 2. The situation is very 
different in the case of type T5 spectrum: near the critical 
density, Sform is insensitive to J, and the critical value of 
J21 is determined almost entirely by H2 photo-dissociation. 

As mentioned above, our J21'' values are a factor of 
3-10 lower than found by OMOl and other previous stud- 
ies. The one-zone models can also be used to identify the 
reason for this discrepancy. In particular, we have system- 
atically varied, one-by-one, each of the five important reac- 
tion rates enumerated above, and re-computed J2"* as these 
rates were varied. We found that the difference in Jfi'' is al- 
most fully accounted for by the difference in our adopted 
H2 colhsional dissociation rate (fcis). In particular, when we 
change our fcis to be the same as in the one-zone calcula- 
tions of OMOl, but leave our other rates unchanged, we find 
Jir* = 1.1 X 10^ (for the type T4 background) and 3.4 x 10* 
(for the type T5 background). For the T4 case, this new 
J2"* is very close to that found by OMOl. In the T5 case, 
our J2"* is still a factor of ~ 3 lower. We have varied all 
of the other important reaction rates (i.e., setting them to 
the values used in OMOl), but found that this discrepancy 
remains. We suspect the remaining difference can be due to 
the different cooling functions used. 

In the Enzo runs, as well as in our one-zone models, 
we adopted the H2 colhsional dissociation rate from Martin 
et al. (1996), while OMOl used a rate based on papers by 
Lepp & ShuU (1983), Shapiro & Kang (1987) and Palla et 
al. (1983). In Figure[8l we show fci5 used in OMOl and Enzo, 
as a function of temperature, for H densities of 10~^ cm"'' 
(low density), 10'' cm""' and 10^ cm~''(high density). As 
this figure shows, the rate used in Enzo is higher than the 
rate used in OMOl, especially at the intermediate density, 
which is close to the critical density, and where the differ- 
ence is more than a factor of 10. This large difference at 
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Figure 8. The H2 coUisional dissociation rate (fcis) adopted by 
Enzo and OMOl, as a function of temperature, at three differ- 
ent densities. Thick curves: rates from Martin et al. (1986), used 
by Enzo; Thin curves: rates used by OMOl. Solid, dashed and 
dotted curves show the rates at nj] = 10"'^, 10'^ and 10^ cm~'^, 
respectively. 



temperature to 6000 K at density lower than ~ 3x 10^ cm~"^. 
This is consistent with Figures [3] and |4l where we see that 
the temperature drops at larger radius (lower density) in 
the case of J21 — 10^ and the T4 spectrum than in the case 
of J21 = 3 X 10^ and the T5 spectrum. Figures |6] and [7] 
also show that the minimum in the total dissociation rate 
occurs at a higher density in the case of the T5 spectrum. 
As Figure [8] shows, kis decreases significantly - by about 
one order of magnitude - when the temperature is lowered 
to from 8000 K to 6000 K. 

The H2 formation rate (fcio) and the H~ formation rate 
(kg) also dependent on temperature, but they work in the 
wrong direction to explain the increased J^"* for halo C: 
fcg and kio axe lower when the temperature is lower, which 
would tend to reduce J21''. However, the temperature depen- 
dence of fcg and fcio are both much weaker than that of fcis. 
When the temperature is changed from 8000 K to 6000 K, fcg 
and fcio decrease by 23% and 5%, respectively. We conclude 
that random temperature variations at the level of ~ 20%, 
which naturally occur due to variations in the strengths of 
shocks that occur in the collapsing gas, can account for the 
scatter we observed in J"^*'. 



the intermediate density is primarily due to the difference 
in Tier. OMOl adopted rid from Lepp & ShuU (1983), which 
considered only vibrational transitions. Martin et al. (1996) 
takes into account all rotational and vibrational-rotational 
transitions, and suggest that accounting for both types of 
transitions reduces ricv by a factor of ~ 10. This then ex- 
plains and justifies the reason that our J21'* values are lower 
than in OMOl. 



3.4 The Origin of Scatter in J2T' 

Another interesting result from the Enzo runs is that J21'* 
found in halos B and C are higher than both in halo A and 
in the one-zone model. This might be due to variations in 
the gas temperatures in the collapsing halos. During the as- 
sembly of the halos, due to the lack of spherical symmetry 
and variations in the merger histories, shocks occur at var- 
ious densities and Mach numbers. This is known to cause 
a ~ 20% scatter in the temperature at fixed density (e.g. 
Loken et al. 2002). The comparison of Figures [T] and [5] with 
Figures [3] and |4] shows that the temperature at rid in halo 
A is always higher than in halo C for the J2i's at which 
these two halos differ (J21 = 10^ for type T4 background, 
J21 = 3 X 10* for type T5 background). Among the five most 
important reactions discussed above, H2 coUisional dissoci- 
ation is most sensitive to temperature. If the temperature 
is lower, the coUisional dissociation rate is also lower, which 
requires a higher J21'* to compensate. 

We used the one-zone model (which, recall, gives good 
agreement with halo A) to check that this is the main reason 
for the different J*^"* values for halos A and C. We find that 
J|;'* indeed becomes slightly larger than 10^ (3.2 x 10'') 
in the case of a T4 (T5) background if the temperature is 
artificially decreased by « 20%, from f=s 8, 000 K to 6, 000 K 
at a density of 5 x 10^ cm"^ (5 x 10^ cm"^). Note that the 
temperature difference occurs at a higher density in the case 
of a T5 background. Jf"' could not reach 3 x 10* if we set the 



3.5 H2 Self-Shielding 

The value of J21'' is directly related to the self-shielding 
factor /sh that we computed with an approximate method 
(equations |4] and [Sj . We here check the accuracy of these 
approximations a-posteriori, using the simulation outputs. 
We computed A'h2 and fsh by integrating the H2 profile 
from the outside in, and compared these values with those 
obtained from Equations |4] and [S] 

We found results that were qualitatively similar in all 
the cases, and show only one example of such a comparison 
in Figure |9l for Halo A in the presence of a T5 background. 
The thick (thin) curves show the results from the numerical 
integration (approximate method). Since the halos are not 
spherically symmetric, integrating along different sight-lines 
gives different results. In our numerical integration, we first 
obtain a spherically average halo profile from the simula- 
tions outputs. At a given density n, we then obtain /ah, by 
averaging the self-shielding factor, over all directions, at the 
radius with this density. For A^H2 , simple averaging makes 
little sense (the photodissociation rate is linearly propor- 
tional to fsh, but not to A''h2), so we only show the results 
of integrating along the single radial sight-line in the direc- 
tion away from the halo center. 

Figure [9] shows that over the interesting density range, 
7Vh2 and fsh obtained from Equations |4] and [5] agree within 
a factor ~ 10, with the values obtained from the non-local 
integrations. In those cases when the halo gas cools, A''h2 
computed with the approximate method is always larger, 
by a factor of a few, while in those cases when the halo 
stays hot, A^H2 computed with the approximate method is 
larger at low density, but becomes very close to the non- 
local value at densities above ^ 10^ cm~^. This is explained 
simply by the fact that Equation [5] assumes that the H2 
fraction is constant. As Figures [T] - 14] show, this is a poor 
assumption for "hot" halos and in the outer region of "cool" 
halos, where /h2 decreases with radius. Only in the inner 
region (R <.10 pc) of "cool" halos does /h2 become approx- 
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Figure 9. The H2 column density and self-sliielding factor in the 
Lyman— Werner bands for Halo A and a T5 background. The thin 
curves were computed from the approximations in Equations |4] 
and |5] the thick curves were obtained by numerically integrating 
the H2 profile in the simulations. The various line types corre- 
spond to different intensities J, as labeled. See text for details. 



imately constant (~ 10^**). From any given location, the 
sightline directly pointing outward (away from the halo cen- 
ter) is typically the least self-shielded, and the mean non- 
local self-shielding factor is therefore lower than one would 
obtain assuming a spatially constant (except close to the 
halo center). As a result, the thick curves in the /ah panel 
of Figure |9] are offset downward from the positions expected 
from the panel using Equation |4] 

Since at Ucr, /h2 computed with the approximate 
method is always slightly lower than its non-local value, we 
expect that J|"' would be slightly lower if /hj were com- 
puted exactly. Furthermore, self-shielding could be overesti- 
mated in our treatment, since we have ignored bulk motions 
of the gas. In the extreme limit where we set the H2 self- 
shielding factors to zero, in the one-zone model, we find 
that the critical fluxes are reduced by factors of « 2 and 
« 30, to J^r* = 19 and J^j" = 440, for the type T4 and T5 
backgrounds, respectively. 



4 THE MASS OF THE CENTRAL OBJECT 

The thermodynamical properties of the collapsing gas have 
a large impact on the final object that forms in the core 
of the halo. An estimate for the mass of the central ob- 
ject, in particular, can be obtained by the following simple 
argument, under the assumption that the object is a sin- 
gle (super-) massive star. There exists a radius at which the 
mass accretion time-scale face equals the Kelvin-Helmholtz 
time scale txH for a proto-star, with the proto-stellar mass 
equal to the gas mass enclosed within this radius. The gas 
outside this radius does not have the time to be incorpo- 
rated onto the proto-star before it settles to the zero-age 
main sequence (ZAMS); conversely, the mass inside this ra- 
dius will accrete onto the protostar before the star settles 
to the ZAMS (e.g. Abel et al. 2002; Omukai & Palla 2003; 
O'Shea & Norman 2007). While this argument ignores var- 
ious feedback processes that can occur and affect the final 
stellar mass, it gives a useful order-of-magnitude estimate. 
See McKee & Tan (2008) for a more detailed discussion in 
the context of first star formation in minihalos. 

In Figures fTUl and [TTI we show the accretion time scale 
tacc = R/{vr) as a function of the gas mass Afgas enclosed 
within a sphere of radius R. Here (vr) is the mean radial 
in-fall velocity at R. These figures show that at enclosed 
masses of < 10^ M©, the accretion rates in the "hot" halos 
are 1-2 orders of magnitude higher than for the "cool" halos. 
The Kelvin-Helmholtz contraction time scale for metal-free 
stars is approximately 10^ years, with a relatively weak de- 
pendence on the proto-stellar mass (Schaerer 2002). Within 
this time, the central object in a "hot" halo will accumulate 
10^ Mq , while the object in a "cool" halo could accumulates 
lO'^ — 10** Mq. The latter values are very similar to previ- 
ous estimates for the masses of the first stars that form in 
minihalos (Abel et al. 2002; Bromm et al. 2002; O'Shea & 
Norman 2007); this is not surprising, given that the densities 
and temperatures reached in these "cool" halos via efficient 
H2-cooling are very similar to those in lower-mass miniha- 
los. However, the value 10^ Mq for the "hot" halos is much 
larger; in particular, a supermassive, metal-free star with 
this large a mass suffers from post-Newtonian instabilities, 
and ultimately collapses without an explosion to produce a 
SMBH of the same mass (Fuller, Woosley & Weaver 1986). 

Given that the effect of the UV background is to dis- 
sociate molecules, and to reduce the efficiency of cooling, it 
may be surprising that the gas accretion rate in the "hot" 
halos is higher than in the "cool" halos, whose gas cools 
more efficiently. However, we note that the accretion rates 
in Figure [10] and [11] are shown at different times: although 
the gas collapse in the "hot" halos is more rapid, it occurs af- 
ter a significant delay. Once the collapse begins, it proceeds 
over the dynamical time scale, 

Here, G is the gravitational constant, and p is the total (gas 
+ dark matter) density. On the other hand, the infall speed 
is modulated by the sound speed. Once the infall becomes 
supersonic, weak shocks occur, which tend to slow the infall. 
These shocks also provide an additional source of heating 
the gas, thereby elevating the sound speed. As a result, the 
sound speed and the infall speed tend to trace each other. 
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In Figure 1121 we show spherically averaged profiles of the 
radial velocity (solid curves) and the sound speed (dashed 
curves) in halo A (and type T4 backgrounds) . The thick and 
thin curves are for J21 = 10^ and J21 = 10'', respectively. 
Although not exactly equal, the sound speed and the radial 
infall speed are always of the same. A similar conclusion 
could be drawn from Figure 1131 where we show three time 
scales: the cooling time (solid curves), the sound crossing 
time (defined as R/cs \ dashed curves) and the free-fall time 
(dotted curves). The thick and thin curves are for the same 
two values of J21 as in Figure 1121 Inside the virial radius 
(i.e., within a few pc, where the gas is first shocked and 
heated to 10*K) the cooling and sound crossing times both 
trace the free-fall time scale to within a factor of ~ 2. It is 
particularly revealing that the good match between the three 
time-scales holds over a factor of ~ 300 in radii. Using the 
condition that the sound crossing time equals the free-fall 
time, 



R/ Csound — R 



32Gp' 



(16) 



together with M 
and T, we find 



pR^ to convert p in Equation [15] to M 



(17) 



where a numerical factor of order unity (which depends on 
the actual density profile) has been omitted. Note that to 
within these factors, this result agrees with the well-known 
result for the accretion time-scale in the self-similar collapse 
of a singular isothermal sphere (Shu 1977). In their simula- 
tions, O'Shea & Norman (2007) find the same scaling of the 
accretion rate with sound-speed, and attribute this scaling 
to the physics of the collapsing isothermal sphere; note that 
our explanation - namely that weak shocks limit the infall 
velocity to be close to the sound speed - is somewhat differ- 
ent. Equation [T7] predicts a linear dependence between taccr 
and the enclosed mass A/, as well as tacc oc y~3/2. jj^^j^ 
of these scalings are in good agreement with the accretion 
times shown in Figures [TOl and [TTI As discussed in §[3]above, 
the central gas temperatures in the "hot" halos are a fac- 
tor of ~ 6 higher than those in the "cool' halos (with the 
exception of the runs with J21 ^ 1). Based on the above 
argument, tacc is therefore expected to be a factor of ~ 15 
smaller, at a given Mgas, for the "hot" halos. This factor is 
indeed in good agreement with the difference in the accre- 
tion rates for the two types of halos shown in Figures [TOI and 

El 



5 CONCLUSIONS 

By performing a series of simulations with the AMR code 
Enzo, we determined the critical intensity of the UV back- 
ground needed to suppress the H2 cooling in dark matter 
halos with virial temperatures of Tvir >10^K. We have mod- 
ified Enzo to include H~ photo-dissociation, as well as self- 
shielding of H2 in the Lyman-Werner bands. To increase 
the computing speed, the self-shielding effect was calculated 
with an approximate method that uses only local physical 
quantities; we have shown, however, that this simple approx- 
imation is accurate to better than an order of magnitude. 




J2i=1 X 10" 
J2i=1 X 10 
J2i=3x10; 
J2i=1 xlO^ 
J21=1 X 10-* 



10° 



Figure 10. The local accretion time-scale tacc as a function of 
the enclosed gas mass Afgas for Halo A and a T4 background with 
different intensities, as labeled. 




10 IC 10° 

Mgas (Me) 

Figure 11. Same as Figure [TOl except for type T5 backgrounds. 



10° 




Figure 12. The radial velocity (solid curves) and sound speed 
( dashed curves) as a function of radius in halo A in the presence of 
type T4 backgrounds. The thick and thin curves are for J21 = 10^ 
and 10^, respectively. 
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R(pc) 

Figure 13. The cooling {solid curves), sound-crossing (dashed 
curves) and free— fall time scale {dotted curves) as a function of 
radius in halo A, in the presence of type T4 backgrounds. The 
thick and thin curves are for J21 = 10^ and 10^, respectively. 



In agreement with earlier results from a one-zone model 
(OMOl), we found that, depending on whether the back- 
ground intensity is above or below the critical value, the 
halos can be unambiguously grouped into two categories: 
"hot" halos and "cool" halos. In general, "hot" halos have 
much higher central gas temperature, density, and electron 
fraction and a much lower H2 fraction. In the "cool" halos, 
the temperature drops inside 7? ~ 10 pc, which corresponds 
to a local density n ~ 10^ cm"''. 

For type T4 (T5) spectra, J|"* is found to be in the 
range of 3 x 10^ - 10^ (10* - 3 x 10'') for halo A, and 10^ 
- 2 X 10^ (3 X 10* - 10^) for halos B and C. These values 
show that the critical flux varies significantly from halo to 
halo; we attribute this variation to scatter in temperature at 
a given density. Most importantly, the values of J21'* we find 
are a factor of 3-10 lower than previously estimated. Using 
one-zone models with simplified dynamics, we studied the 
dependence of J21'' on each of the important reaction rates. 
We have shown that the difference between our results and 
previous estimates can be attributed to our adoption of a 
different, more accurate, H2 coUisional dissociation rate. 

In those halos in which H2 cooling is suppressed, the 
gas cools efficiently, but remains relatively hot, at a tem- 
perature near T 8000 K. While gas collapse starts with 
a significant delay, we have shown that, as a result of the 
elevated temperature, the gas accretion rate in these halos 
is ultimately increased by « 1 — 2 orders of magnitude. As 
a result, a supermassive star with a mass of 10^ Mq may 
form in the cores of these halos (compared to 10'^"'' M© 
stars in the presence of H2 cooling), ultimately producing a 
supermassive black hole (SMBH) with a comparable mass. 

The critical UV fiuxes required to suppress H2 cooling 
m halos withTvir >10*K are high compared to the expected 
level of the cosmic UV background at high redshifts. Most 
likely, the halos exposed to a super-critical UV ffux are a 
small subset of all Tvir ^IC'K halos that happen to sam- 
ple the bright-end tail of the fluctuating cosmic UV back- 
ground. This makes the reductions in J2"* that we have 
found particularly significant: according to a model for the 



fluctuating UV background, as sampled by DM halos (Di- 
jkstra et al. 2008) the 10-fold decrease in J"^^ can increase 
the number of candidate DM halos, where direct SMBH for- 
mation may be feasible, by a factor of « 10^. 
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Reaction 




Rate CoefBcient k (cm^ s ^) or Cross-section cr (cm^) 


(1) H + e- H+ + 2e- 


fei = exp(-32.71396786375 






+13.53655609057(lnTev) 






-5.739328757388(lnTcv)^ 






-|-1.563154982022(lnTcv)^ 






-0.2877056004391 (InTev 






-|-0.03482559773736999{lnT'cv)'^ 






-0.00263197617559(lnTcv)'^ 






-1-0 0001 1 lQ5439538fil HnT ir"!^ 






-2.039149852002 x 10-«(lnTev)*) 


(2) He + e" 


He+ + 2e- 


fe2 = exp(-44.09864886561001 

-|-23.91596563469(lnTev) 

-10.75323019821(lnTev)^ 

-|-3.058038757198(lnTcv)^ 

-0.5685118909884001(lnrov)'' 

-l-0.06795391233790001{lnT<,v)^ 

-0.005009056101857001(lnTev)'^ 

-1-0 00020fi72,36157,507('lnT -s/"!^ 

-3.649161410833 x 10-6(lnTev)®) 


(3) He+ + e 


- -* He++ + 2e- 


fes = exp(-68.71040990212001 

+43.93347632635(lnTev) 

-18.48066993568(lnTev)^ 

-|-4.701626486759002(lnTev)^ 

-0.7692466334492(lnTcv)*' 

-|-0.08113042097303(lnTcv)^ 

-0.005324020628287001(lnrev)^ 

-|-0.0001975705312221(InTev)'^ 

— 3 Ifi55810fi5fifi5 X 1 0~^C1tiT Ar'l®'! 


(4) H+ + e" 


-» H + hi/ 


ki = exp(-28.61303380689232 
-0.7241125657826851(lnTev) 
-0.02026044731984691(lnTcv)^ 
-0.002380861877349834(lnTev)^ 

-0.0003212605213188796(lnrev)'' 
-0.00001421502914054107(lnTcv)-''' 
-1-4.989108920299513 x 10"''(lnTov)'' 
-1-5.755614137575758 x lO-^'OnTev)'^ 
-1.856767039775261 x lO-^OnTev)* 
-3.071135243196595 x 10-3(lnTev)^) 


(5) He+ + e 


- -» He + hi/ 


fes = 1.54 X 10-9(1 -|-0.3/exp(8.099328789667/Tev)) 
/(exp(40.49664394833662/rev) x T^y^) 
+3.92 x 10-i3/r?y6353 


(6) He++ + 


Hc+ + hu 


ke = 3.36 x 10-i"T-5 (T/1000)-0-2(l + (T x lO-^)"-'')-! 


(7) H + H+ 


^ H^ + ht/ 


kr = 1.85 X 10-23 X T^-^ T ^ 6.7 x 10^ K 




kr = 5.81 x io-i6(T/56200)(-0-6«57iog(T/6620o)) T > 6.7 x 10^ K 



APPENDIX: REACTION RATES AND CROSS SECTIONS 



T and T^y are the gas temperature in units of K and eV, respectively. 
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Reaction 



Rate CoefBcient k (cm^ s ^) or Cross-section cr (cm^) 



H 



(8) H+ 

(9) H + e- 

(10) H + H 



> H2 + H+ 
H- + hi^ 
■^Ha+e- 



(11) H++e- - 

(12) H+ + H- 

(13) H- + H+ 

(14) H2 +e- - 

(15) Ha + H^ 

(16) H2+H2 - 

(17) H2 + H+ - 



■ 2H 



-»H2 + H 

2H 
> H + H- 
3H 

» H2 + 2H 
-»H++H 



(18) H2 - 

(19) H- 



2H- 



2e- 



(20) H- + H ^ 2H + e" 



(21) H- +H+ 

(22) H + hu-^ 

(23) He + hi/ - 

(24) He+ + hu 

(25) H- + hi/ - 

(26) H J + hi/ - 

(27) H2 + hi/ - 

(28) H2 + hi/ - 



■H++e- 

H+ +e- 

* He+ + e- 
Hc++ -f 
-+H + C- 
H-I-H+ 
.H++e- 
•> 2H 



6.0 X 10"^" 

fcg = 6.775 X W-^-'T^y^™ 
kio = 1.43 X 10-9 
fcio = cxp(-20.06913897587003 
+0.2289800603272916(hiTev) 
+0.03599837721023835(lnTcv)^ 
-0.004555120027032095(hircv)-^ 
-0.00031051 15447124016(lnrev)* 
+0.0001073294010367247(lnTev)^ 
-8.36671960467864 x 10-«(lnTev)*' 
+2.238306228891639 x lO-'^(lnTev)^) 
fcii = 1.0 X 10-* 
fell = 1.32 X io-eT-° '^6 

ki2 = 5.00 X 

kvi = 6.5 X lO-^T 
fcl4 = 0.0 (not used in Enzo) 
see expression in Martin ct al. (1996) 
fcie = 0.0 (not used in Enzo) 
kn = exp(-24.24914687731536 
+3.400824447095291 (InTev ) 
-3.898003964650152(lnTev)^ 
+2.045587822403071(InTev)^ 
-0.5416182856220388(lnrev)'' 
+0.0841077503763412(lnTcv)'^ 
-0.007879026154483455(lnTcv)'' 
+0.0004138398421504563(lnTov)'^ 
-9.36345888928611 x 10-''(lnTev)®) 
fcl8 = 5.6 X 10-"cxp(-102124/T)r°-S 
ki9 = exp(-18.01849334273 
+2.360852208681(lnTev) 
-0.2827443061704(lnTcv)^ 
+0.01623316639567(lnTov)^ 
-0.03365012031362999(lnrev)* 
+0.01178329782711(lnTov)'^ 
-0.001656194699504(lnT'cv)'^ 
+0.0001068275202678(lnrev)'^ 
-2.631285809207 x lO-S(lnrev)*) 

k20 = 2.56 X 10-9T^1/8186 

feao = exp(-20.37260896533324 
+1 . 139449335841631 (InTev ) 
-0.1421013521554148(lnrev)^ 

+0.00846445538663(lnTev)^ 
-0.0014327641212992(lnTev)* 
+0.0002012250284791(lnrev)'''' 
+0.0000866396324309(lnrcv)*' 
-0.00002585009680264(lnTev)^ 
+2.4555011970392 x 10-'5(lnrev)* 
-8.06838246118 x 10-*(lnTev)^) 
fe2i = 4 x 10-''T-i-4exp(-15100/T) 
fe2i = 1 x 10-8T-0-4 

<T22 = 6.3 X 10-l«(i^j^)4exp(4 - 4(tan-l)e/e)/[l - 
"■23 = 

(724 = 1.575 X 10 

k2B = 10-"'tt,/2i; a = 2000 for T4 spectrum, 0.1 for T5 spectrum 
see expression in Shapiro & Kang (1987) 
see expression in Shapiro &; Kang (1987) 

feas = 10-12/3 J21; /3 = 3 for T4 spectrum, 0.9 for T5 spectrum 



Tev ^ 0.1 

TeV > 0.1 



T ^ 617 K 
T > 617 K 



TeV < 0.1 
TeV > 0.1 



T < 10* K 
T > lO'' K 



0.694 X 10-i8[(^)i-S2 + (^)3. 23^-1 
54.4cV 



exp(-27r/e)];6 = ^j|i^-l 



-IS 



(54^)4 ,^p(4 _ 4(tan-l)e/e)/[l - exp(-27r/£)]; 6 = / 



hu 
54.4 



